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Dedicated to J. Zittartz on the occasion of his 60th birthday 

Abstract. In order to study the effect of interaction and lattice distortion on quantum coherence 
in one- dimensional Fermi systems, we calculate the ground state energy and the phase sensitivity 
of a ring of interacting spinless fermions on a dimerized lattice. Our numerical DMRG studies, in 
which we keep up to 1000 states for systems of about 100 sites, are supplemented by analytical con- 
siderations using bosonization techniques. We find a delocalized phase for an attractive interaction, 
which differs from that obtained for random lattice distortions. The extension of this delocalized 
phase depends strongly on the dimerization induced modification of the interaction. Taking into 
account the harmonic lattice energy, we find a dimerized ground state for a repulsive interaction 
only. The dimerization is suppressed at half filling, when the correlation gap becomes large. 

PACS. 71.10.-w Theories and models of many electron systems - 75.10.Jm Quantized spin models 



1 Introduction 

Recent experiments on CuGeC>3 [|J and NaV2C>5 J|,||, 
which show at low temperature a transition to a non- 
magnetic ground state accompanied by a structural tran- 
sition, led to renewed interest in spin-Peierls systems, as 
well as in the more general question of structure ver- 
sus correlation induced metal- insulator transitions. The 
model of spinless fermions, which we consider here in de- 
tail, describes certain aspects of both, and it contains, as a 
special point in parameter space, the experimentally rel- 
evant isotropic Heisenberg antiferromagnet. In addition, 
■ the model contains, in the limit of an Ising type interac- 
tion, a szenario which is closer to that found for a more 
realistic model, including the electron spin, namely the 
Hubbard-Peierls model. Here we concentrate on the spin- 
less fermion model, which nevertheless shows a surpris- 
ingly rich phase diagram (at zero temperature) as a func- 
tion of interaction and dimerization. 

In the next section, we introduce the various repre- 
sentations of the model. In Sect. 3, we quantitatively de- 
termine the extension of the delocalized phase, which was 
predicted for an intermediate attractive interaction. There- 
by we identify the region where the fermion-phonon cou- 
pling is irrelevant, i.e. the ground state remains extended. 
In Sect. 4, we consider a repulsive interaction, where we 
find that a stable dimerization develops. For a very strong 
interaction the dimerization is reduced again, because con- 
flicting ordering occurs. In the summary, Sect. 5, we present 
the u-V phase diagram, as well as the results for the dimer- 
ized state. 



2 The model 



As a starting point for the study of a general spin-Peierls 
system in one dimension, we consider an XXZ model with 
a dimerized interaction in a local magnetic field; the latter 
can be considered as a staggered magnetic field resulting 
from the surrounding chains H , or resulting from random 
magnetic impurities. In the following, we concentrate on 
the dimerization-induced modification of the interaction, 
and, for comparison, include the magnetic field only occa- 
sionally. Furthermore, we neglect frustration effects which 
arise from next-nearest-neighbor couplings. Thus we start 
with the following Hamiltonian: 



M 



A I 



71=1 



Ko 2 



(1) 



where J n (u) = J(l — (— ) n Xju); the local lattice distor- 
tion, u n , is assumed to be given by u n — 2(x n — %n)/ a = 
(-l) n u. For the "clean" XXZ model, i.e. for u = Q, h n = 0, 
one finds three phases: a ferromagnetic phase for A > 1, a 
gapless phase for — 1 < A < 1, whose low lying excitations 
are given by those of a Luttinger liquid, and an antiferro- 
magnetic phase for A < — 1. The spin model is experimen- 
tally most relevant for A — — 1, i.e. the isotropic antifer- 
romagnetic Heisenberg case. The corresponding fermion 
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model is obtained via the Jordan Wigner transformation: 3 Phase transition at attractive interaction 



2n„ 



1, 



It is customary to change the notation: J = t, J A = 
— V/2, and e m — —2h m '-> furthermore we now assume that 
the coupling to (static) phonons can be varied indepen- 
dently in the hopping and the interaction. The result is 
the following: 

-fffcrmion = — ^ ] t m (u) (Cm c m+1 + c m+l c rn) 
m 

m 

+ ^2e m n m + M^-u 2 , (3) 

m 

where t m (u) = t(l-(-) m \ t u), V m (u) = V(l-(-) m \ v u). 
Furthermore, we consider twisted boundary conditions, 
Co = e^cm- The length of the chain is L — Ma, the num- 
ber of electrons N, and the filling no = N/M. In addition, 
we restrict ourselves to half filling, and set t = 1 and 
A t = 1 in some of the formulas below. The corresponding 
instabilities in the fermion model are at V = — 2, towards 
a phase separated state, and at V — 2, to a charge den- 
sity wave (CDW) state. (Note that at V = 2, the Ak F - 
backscattering process becomes relevant [|).) Numerically 
this CDW state is difficult to determine, because, for finite 
systems, the hopping lifts the twofold degeneracy of the 
ground state via a symmetric superposition, and a uniform 
density is found. For the analytic considerations, we use 
the fact that the system is a Luttinger liquid in the gap- 
less phase, and that the dimerization can be considered a 
perturbation. In the bosonized form ||, the Hamiltonian 
can be written as follows: 

p { V -[ dxi p{ x ,t)f + -[d tV >(x,t)f 
2tt l_ g vg 

2tt\u . . . .. nVa ,. , ..1 

■ sin [2cp(x, t)] - — cos [4tp(x, t)] } ,(4) 



where A = A t £ + XyVno/2, and Iq is a short distance 
cut-off, w 2a. Only the first order correction to the 2kp- 
process, oc Vu, has been included. The Fermi velocity is 
given by v — [nt sin(2n)] / (tt — 2i]), and g — tt/^t], where rj 
parameterizes the interaction according to V= —2t cos(2?y). 
Furthermore, the density is given by 

\<=x/a = Pi x ) - n a = — + j- cos (2k F x + 2ip). (5) 

Z ' IT If) 

Clearly, the term oc cos [Aip(x,t)j, describing Akp scatter- 
ing processes, is responsible for the CDW ordering. For 
completeness we like to state that some properties of the 
dimerized Heisenberg chain, especially the spectrum, can 
also be obtained j7| from the equivalent relativistic model, 
the massive Thirring model. 

The numerical results, presented below, are obtained 
on the basis of the XXZ model, Eq. (flh. 



The first conjecture that the dimerized XXZ model under- 
goes a transition from the dimerized spin singlet phase to a 
phase with free spins at A = V2/2, was given in || based 
on a mapping onto the Ashkin- Teller-model. There, an es- 
timate of the phase boundary, based on a series analysis, 
has been given. (Note the different sign convention.) For 
the XXZ model in a staggered magnetic field, Alcaraz || 
again found this critical point on the basis of predictions of 
conformal field theory. Thus it seems that the periodic dis- 
tortion of the <7 2 c7 2 -interaction does not change the point 
of the transition. In the case of the bosonic model, assum- 
ing an interaction of the form U n cos nip, one finds the 
following renormalization group equations [fill : 



— = U j2--g 



E 

ni+ri2— n 



U ni U n2 . (6) 



The critical value or a model containing only a single pe- 
riodic term, g = 8/n 2 , is in accordance with the exact 
solution of the sine-Gordon model Jll],[[D ■ One concludes 
that for the "clean" case (u = 0), there exists an instabil- 
ity for g = 1/2, i.e. V = 2. Consequently, considering the 
sin(2<^)-contribution in (0), we expect interesting effects 
for g — 2, i.e. V — — v2, provided the dimerization, u, 
assumes a finite value. 

Applied to our model, the critical behavior near the 
transition between a localized and a delocalized state is 
determined by the Berezinskii-Kostcrlitz-Thouless equa- 
tions ill: 



du 
~dl 



u(2 - g) 



dg 
dl 



2 \ 2 2 

-g A u 



(7) 

(8) 



Considering the bosonized Hamiltonian, the equivalence 
of the "staggered" and "dimerized" XXZ model is also 
apparent: For the dimerized model, the relevant operator 
is proportional to u sin 2ip, while for a staggered field, it is 
cx h cos 2ip. The phase shift of 7r/2 corresponds to the cou- 
pling to the bond (in the hopping term) versus coupling 
to the site. As a consequence, the ground state for the 
dimerized model is a spin singlet, with ip m - m — ir/2 and 
ofj = 0, while in a staggered magnetic field, an antiferro- 
magnetic ground state, with ip min = and a 2 = (—1)™, is 
obtained. The critical behavior, however, is the same for 
both models. 

We numerically verify these predictions and also de- 
termine quantitatively the phase boundary between the 
localized and the delocalized phase. In addition to the 
critical point at u = 0, g = 2, physical arguments suggest 
that the delocalized phase should not exceed the point 
u = 1, V = — 2: the transition at V — —2, being first 
order, is stable against perturbations ||, and for \ t u = 1 
every second bond is cut and therefore the system clearly 
is localized. Compared with the random disorder case [jlij , 
one is tempted to expect - incorrectly - that the delocal- 
ized phase is smaller. For the random case, this is due to a 
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stronger renormalization of g, namely dg/dl — — 0.5g 2 W 
jL{|; compare this with equation @. 

We will use the phase sensitivity, i.e. the reaction of 
the system to a change in the boundary condition, to de- 
termine this transition numerically for systems with fi- 
nite size. The ground state energy E{<j>) depends on the 
boundary condition, expressed via the phase <fi. In partic- 
ular, we determine below the energy difference between 
periodic and anti-periodic boundary conditions, AE — 
(—) N (E(0) — E(n)). Recall that, for a clean system, the 
ground state energy and the finite size corrections can 
be obtained from the Bethe Ansatz p^] . At half filling 
(and for odd particle number), the result, in the Luttinger 
regime, is given by Q 



E M (4>) - MSoa = 



TTV 

6M 



l-ig- 



(9) 



where is the energy density in the thermodynamic 
limit. Thus MAE = irvg/2, independent of M, for the 
metallic state. In an insulator, on the other hand, the sys- 
tem cannot react to a twist in the boundary condition, i.e. 
MAE is expected to decrease with system size. 

For non-interacting fermions, the finite size corrections 
for the Su-Schrieffer-Heger (SSH) model jig] , (u) oc 
u 2 lnu, can be determined in an elementary way, with the 
help of the Euler-McLaurin formula. However, it is neces- 
sary to handle the van Hove singularity at the band edge 
carefully for periodic boundary conditions, while for anti- 
periodic boundary conditions, this term is absent. With 
h = 2tt/M, we obtain the following result: 



E M (u) - M £oo (u) = ( — lnu - — ln(V« 2 + >i 2 - /l ) 



1 1 h 2 1 -u 2 ) 1 , 
-u-\ , h 

2 12 ^h 2 + u 2 12 



(10) 



Thus, for small systems, h 3> u, we obtain a correction 
linear in u which does not depend on the system size, and 
a correction to the u 2 In it-term: 

1 u 2 j u 

which we plot in figure |l| For large systems, h -C u, the 
corrections are of order 1/M: 



12 



12 u ' 



Thus, for small enough systems, or small enough u, the 
phase sensitivity is MAE(u) oc —uM; this behavior re- 
sults from the zone boundary contribution. 



For V 7^ 0, we combine the treatment of Loss 19 of a 
Luttinger liquid on a closed ring with the scaling equation 
(Q), and a first order perturbative calculation, with the 
following result: 



MAE(Xu) = MAE(Xu = 0) - Xu 



2vr / M 



v \M 



2-9 



(11) 



This is similar to the case of a single impurity jwj, but 
valid in a restricted regime only. This estimation cannot 
be extended to large M or it, because, in a finite sys- 
tem with a gap A(u), the "free motion" is only seen if 
L/£ oc MA(u) oc Mu is small. For larger systems, the 
phase sensitivity depends exponentially on the correlation 
length, as discussed in ^0|. In addition, we have neglected 
the renormalization of the interaction parameter g, see 
equation (|J), which is reasonable only for V < 0, as can 
be seen easily by integrating the RG equations (^) and (||) . 
For example, only the samples with u — 0.01, for system 
sizes M — 24 and M = 50, can be fitted with equation 
(|i"T|). Already for u = 0.05 and M — 24 deviations from 
the first order treatment are obvious; compare figure ^. 

Our aim is to determine quantitatively the phase bound- 
ary for the dimerized XXZ model (Ay = At), and to com- 
pare the results with those in the case of Ay = 0. In the 
latter case, one may compare with results obtained from 
the massive Thirring model, or the RG equations. In our 
density matrix renormalization group calculations [ pl[ , we 
kept 1000 states and performed 9 finite lattice sweeps. 
First, consider Ay = 0. In order to determine the phase 
transition near u = 0, we plot in figure |^ the phase sen- 
sitivity for M = 50 and 100, and small u, and confirm 
that the transition occurs near V f» — 1.4. (The transition 
point is defined as the point where the phase sensitivity 
becomes independent of the system size.) Proceeding to 
higher values of u, we can follow the phase boundary, as 
shown in figure [|. Already at u = 0.4, we find a completely 
localized phase, and we verify the presumption that the 
delocalized phase extends to the line V = —2. In figures 
[| and [|, we show the difference between the two cases, 
Ay = A t versus Ay = 0. There is no drastic change in the 
qualitative behavior for small u, and the transition point 
remains, as expected from (Q) and the renormalization 
group equation (0), at V ~ —y/2. We also find, in the case 
Ay = At , that the phase sensitivity decreases more slowly 
for an attractive interaction, and faster for a repulsive in- 
teraction. Increasing the dimerization further, we see that 
the delocalized phase is considerably larger for Ay = At, 
and indeed extends to the point u = 1. In comparison 
with the disordered model p~ 
a stronger interaction, and t 
to the line V = -2. 



the transition occurs at 
re delocalized region extends 



4 Interaction and dimerization 

In the last section, we determined the u-V phase diagram, 
under the assumption that the lattice distortion is a fixed 
parameter. In a more realistic model, the harmonic lattice 
energy will compensate the energy gain of the fermions 
or the spins, respectively, resulting in a finite equilibrium 
value of u. 

We consider the ground state energy for periodic bound- 
ary conditions, studying the lattice effect for — 2 < V < 2 
first. The numerical data (see figures and ||), obtained 
for 100- or 200-site systems, while keeping about 300 to 
400 states (for about 100 sites the finite size corrections 
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are very small and the term linear in u is absent), show 
that we can distinguish three regimes. For V between —2 
and 0, the gain in the fermionic energy grows parabol- 
ically, AE(Xu) oc -(Au) 2 . For < V < 2, the energy 
gain is stronger, oc u a , with an exponent which reaches 
its minimum of 4/3 at V — 2, where correlations change 
the behavior. 

Explicitly, from an analysis of the massive Thirring 
model, one finds Q the energy gain, AE(\u) — E{\u) — 
E(u = 0), to be given by 



c 2 u 2 , 



-V2<V<0: AE(Xu) = c x u a - 
a = 2/(2 - <?) > 2; 
< V < 2 : AE(Xu) = -cm" + c 2 u 2 , 

4/3<ev<2, (12) 

with ci > C2 > 0; compare figure 9. The non-quadratic 
term can also be obtained from a self-consistent treatment 
of the perturbation. In this so-called self-consistent har- 
monic approximation (SCHA) Jll|,^3| an optimal quadratic 
approximation for the periodic potential around one of its 
minima is chosen, which corresponds to the following re- 
placement: 



Ucosntp -> U[l - y (v 2 - (^ 2 ))]e~ n2 <^/ 2 ; 



here, explicitly 



2 Z <3o 



ngU 



A*_ 
v 2 ' 



(13) 



(14) 



The quadratic term in equation (|l^) is related to the tun- 
neling between the minima of the periodic potential . 
The logarithmic dependence on u for V — 0, which was 
found in 18 , is contained in the above formulas: 



V = : AE a = - — (Au) 2 (ln^- - 0.5) 

TT Xu 

MV , 4 

V m : AE = AE (Au) 2 In 2 — . (15) 

7T Xu 

Adding the harmonic lattice energy, MK u 2 /2, we real- 
ize that a finite equilibrium dimerization cannot be be 
stabilized for an attractive interaction. For small |V|, the 
second logarithmic term in (|i~5| ) always leads to a positive 
slope of AE(Xu) at u = 0; and for large (negative) |V|, we 
find, as shown in the previous section, a delocalized phase 
for small u. 

For the non-interacting case (V = 0), the stable dimer- 
ization Xuq and the energy gain are given by well-known 
expressions |M for CDW-systems. With 7 = X 2 t/Ko, the 
results are 



Auo oc e 1 ^ 7 , 
A(Xu ) oc e" 1 ^, 
^totaiCAuo) oc e" 2/7 . 



(16) 



For a repulsive interaction, we find the following algebraic 
dependences: 



. 2 -g 
Alto c< 7 2 - 2 s , 

A(Xuq) oc j^F), 

-Etotal(Auo) oc 7 T= «, 



(17) 



and, for example, 

V = 1 

V = 2 



Aw oc 7 5 ^ 2 , 
Auq oc 7 3 ^ 2 , 



-Etotai(Auo) oc 7 4 ; 

Etotal(Xuo) oc 7 2 . 



Another interesting quantity is the curvature of the to- 
tal energy at the minimum, which is related to the zone 
boundary phonon frequency, L02k F ■ For a decoupled system 
(A = 0), the curvature is given by K , the bare value. The 
fermion-lattice coupling changes this quantity as follows: 



y = 0: E' t ' otal (u )/M = jK 
V>0: E' t ' otal (u )/M = K Q (2-a). 



(18) 



We thus conclude that for interacting fermions, the soften- 
ing of the phonon at the zone boundary does not depend 
on the coupling strength 7, but on the exponent of the 
energy gain of the fermions only. 

Next we consider the regime V > 2, where a correlation 
gap develops, for the "clean" case given by [E5| 



V > 2 : A(V) oc e- 7 * 2 '^ 2 , 
V -> 00 : A(V) oc V. 



(19) 



A gap also develops in a frustrated model, i.e. when next- 
nearest-neighbor coupling is included, for V — 2 and pro- 
vided that Vfrust > Vfrust- This is a consequence of a rele- 
vant 4fc f -backscattering process, and believed to be rele- 
vant for CuGeC>3 [ p6p^ | . As shown in [^8) , a cos 4<p-term 
leads to a Neel (Ising-type) or a dimer (frustration) ground 
state, respectively, depending on its sign. Thus we expect 
that a next-nearest-neighbor interaction will increase the 
dimerization, while an Ising-type interaction reduces it. 
Possible distinctive features between the frustration in- 
duced and the externally, by Xu, driven dimerized state 
are discussed in |29| . 

Considering the numerical data for a strong repulsive 
interaction, V > 2, see figure ^, the energy gain is nearly 
constant at first, when increasing V; only the prefactor 
is slightly reduced. Arguments supporting this observa- 
tion are given below. Increasing the interaction, the en- 
ergy gain is drastically reduced, the exponent approaches 
again 2, and therefore no minimum is found in the to- 
tal energy. To understand this behavior, we consider in 
more detail the extrema of the combined sin-cos-potential, 
U(f) = —VcosAtp — u sin 2<p, where u = —2nvgXu/lo, and 
V — TrvgV/l 2 . For a large interaction, the minima are at 
sin2<po = u/(4V), and the expansion reads 

[%) *-V-£+(w-liy<p- . (20) 
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So for a strong interaction, it is possible to consider the 
M 2 -terms in the potential as a small perturbation to the 
exact results ^5|. Note that, for a strong interaction (i.e. 
the renormalized g is « 0), fluctuations are not important 
since {tp 2 ) « 0. These arguments support the result, ob- 
tained numerically, that the energy gain is quadratic in u. 
For large u, on the other hand, the minima are located at 
cos2</?o = 0, and the expansion for small interaction reads 



U(<p) « (V- - u) + 2{u - W)(<p - Va y 



(21) 



For an interaction V ~ 2, where the effective influence of 
the V cos 4(/3-term is expected to be small, we repeat the 
perturbation analysis, i.e., we evaluate first the it sin 2(/> 



term within the SCHA; thus A(u) 



oc u 



2/3 



Adding then 



the interaction term to the gapped system, and taking into 
account the fluctuations, the dimerization gap is reduced: 



y ' (l + V) ™ 



(22) 



For a frustrating interaction, there is no change in the 
minima, and we find, for a small frustration, a similar 
expression to (p2|), and for a stronger frustration, a linear 
increase in u. 

To check the validity of the perturbation analysis at 
V — > 2, we compare with results from the SCHA and 
from field theory. In the evaluation of the SCHA scheme, 
which considers both perturbations on an equal footing 
following (21), we find the following gap equation, which 
can be solved for g = 1/2: 



A 1 = 



k>A 
2nv 



the result is: 



\2irv J 



A(u) 



4fl 



A(u,V) = — -±= 



(l + V) 2 / 3 ' 



(23) 



(24) 



Clearly, the exponent of the denominator does not coin- 
cide with the one in fl22|). 

However, a more careful treatment, particularly with 
respect to the marginality of cos4c£> at V = 2, introduces 
a logarithmic correction to the dimerization gap Q] , given 

by 



A(u, V) cx 



[l + V In [v/A(u))}V 2 



(25) 



This logarithmic term is a consequence of the scaling 
the marginal operator V(L) — V(l + irV\nL)~ 1 , see 
and therefore is not found within both of the above con- 
sidered approximations. The remaining difference between 
( pi| ) and (|2^), namely the different exponents, can be ex- 
plained by higher order contributions pOfl . Thus the per- 
turbation approach, which led to (|2^), is rather close to the 
correct result, (25), and even the simple SCHA is not too 
far off, as it shows that the interaction essentially changes 
the prefactor, i.e. the gap is ~ A(u). 



The observation that a repulsive interaction in the one- 
dimensional Hubbard model can enhance the dimerization 
|n] can be explained now in the following way: For small 
U , where the charge gap is small, we find an energy gain 
^charge oc -E S pin oc tu 4 ^ 2 ~ 3c ' , with the corresponding in- 
teraction parameters g s = 1 for the spin and g c for the 
charge. At half filling, an intermediate U leads to a large 
enough gap to push down the dimerization in the charge 
channel, but in the spin channel, no gap appears, and the 
system gains energy from the coupling to the lattice. The 
charge gap reduces, nevertheless, the prefactor from t to 
t 2 /U. 



5 Summary 

In summary, we have shown that a periodic distortion of 
the lattice is irrelevant for V < — \/2, resulting in a delo- 
calized phase. The transition point, V c (u — 0) = —V2, 
is not modified by the way of coupling to the lattice. 
However, the extension of the delocalized phase depends 
strongly on the coupling to the interaction, i.e. the param- 
eter Ay. For the spin model (Ay = At), the delocalized 
region is considerably larger than for the case Ay = 0. In 
comparison with a random disorder model jyj, the most 
important difference is the fact that the delocalized region 
extends down to the line V = — 2. The phase diagram is 
summarized in figure [l0|. Due to the difficulties in deter- 
mining precisely the point where the curves for different 
system sizes coincide, we can plot only a rough phase dia- 
gramr; nevertheless the transition seems to be steeper for 
Ay = At, compared to the case Ay = 0. 

Including the harmonic lattice energy, we find that 
the dimerized state is stable, i.e. the total energy devel- 
ops a minimum at a finite uq, for a repulsive interac- 
tion only. With increasing interaction, the dimerization 
increases even further, provided no interaction induced, 
competing ordering occurs. A schematic plot is given in 
figure O. Furthermore, already at V = 4, no equilibrium 
dimerization develops, as can be seen from the prelimi- 
nary data shown in figure [l^, where we plot the gap versus 
inverse system size. (The lines are obtained from a poly- 
nomial fit in X/M for the case u = 0, and an exponential 
fit as mentioned in [fi2|| , respectively.) Obviously, no ad- 
ditional gap from the dimerization is seen for V — 4. In 
contrast, for smaller V, especially V = 2, the dimerization 
contribution is significant. 



We thank Peter Schmitteckert for providing us with the DMRG 
algorithm, and Peter Schwab for useful discussions. This work 
was supported by the Deutsche Forschungsgemeinschaft (For- 
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Schematic plot of the equilibrium dimerization 
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